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We study the zero temperature phase diagram of Ising spin systems in two dimen- 
sions in the presence of competing interactions, long range antiferromagnetic and 
nearest neighbor ferromagnetic of strength J. We first introduce the notion of a 
"corner energy" which shows, when the antiferromagnetic interaction decays faster 
than the fourth power of the distance, that a striped state is favored with respect 
to a checkerboard state when J is close to J c , the transition to the ferromagnetic 
state, i.e., when the length scales of the uniformly magnetized domains become large. 
Next, we perform detailed analytic computations on the energies of the striped and 
checkerboard states in the cases of antiferromagnetic interactions with exponential 
decay and with power law decay r~ p , p > 2, that depend on the Manhattan dis- 
tance instead of the Euclidean distance. We prove that the striped phase is always 
favored compared to the checkerboard phase when the scale of the ground state 
structure is very large. This happens for J < J c if p > 3, and for J sufficiently 
large if 2 < p < 3. Many of our considerations involving rigorous bounds carry 
over to dimensions greater than two and to more general short-range ferromagnetic 
interactions. 
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1. INTRODUCTION 

In this paper we continue our study of the ground state (GS) of lattice spin sys- 
tems with competing ferro (F) and anti-ferro (AF) Ising-like spin interactions. See 
[10-13] for previous results. Such systems are simplified models of real systems 
with both short range attractive interactions and long range dipolar type interac- 
tions. The competitive nature of these interactions is believed to be responsible 
for the formation of mesoscopic periodic structures, such as stripes, in many quasi 
two-dimensional (2D) systems at low temperature, see [1-6, 9, 14, 15, 19-31] for 
several examples of spontaneous pattern formation in physical systems with com- 
peting interactions. See also [16, 32] where such competition is held responsible 
for the development of macroscopic patterns in chemical and biological systems 
described by reaction-diffusion equations. 

While it is simple to understand that the competition between interactions 
acting on different length scales can give rise to mesoscopic structures, it is very 
difficult to predict the optimal shape of these structures. Here we show for a large 
class of interactions that stripes are energetically favorable as compared to other 
natural structures, such as rectangular or square checkerboard. 

The Hamiltonians we consider have the form 

1 x 1 . 

H= -J^ [-J<J| x -y|,i+eu(x-y)](<T x (T y -l) = -^0(x-y)(cr x (7 y -l) , (1.1) 

x^y x^y 

where x e Z. d , cr x = ±1 are Ising spins, J and e are two positive constants (the 
strengths of the F and AF interactions), v is a non-negative potential, symmetric 
with respect to 90° rotations and summable. In the following, we will be mostly 
concerned with d = 2 and v of infinite range. The constants J and e will be 
thought of as being "large" and "small" , respectively. 

The goal is to understand the zero temperature phase diagram as the ratio 
J/e is varied. If e = 0, then the ground state is ferromagnetic. In the opposite 
limit, that is, J = 0, then the ground state displays some non-trivial alternation 
between positively and negatively magnetized spins; e.g., if t>(x) = |x|~ p , p > 
d, then the ground state is the period-2 antiferromagnetic Neel state [8]. As 
the ratio J/e is increased from zero to large values, the GS changes to reduce 
the number of antiferromagnetic bonds, presumably by displaying mesoscopic 
uniformly magnetized structures of larger and larger lengths. It is often assumed 
that the ground state configurations are periodic, and display either checkerboard 
or striped order, depending on the specific choice of the interaction and the 
value of J. In [21], it was shown that for v(x) = |x|~ 3 and J large enough, the 
optimal striped configuration has lower energy than the optimal checkerboard 
configuration. This leads to the conjecture (still unproven) that the ground state 
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configurations of Eq.(l.l) with v (x) = |x|~ 3 and J large display periodic striped 
order. 

There is in fact evidence for the fact that the sequence of transitions to the 
ferromagnetic phase has some universal features [6, 15, 28, 29] and that the emer- 
gence of stripes is essentially independent of the details of the F and AF interac- 
tions. However, the reason for this is still unclear and puzzling because stripes 
break the symmetry of the lattice. 

In this paper, we prove that striped patterns are favored, within a natural class 
of variational states, when the scale of the GS structure is very large compared to 
the range of the FM interaction. A simple explanation of this fact can be based on 
the concept of corner energy, which suggests that the intersection points among 
straight phase separation lines can be thought of as elementary excitations of 
the system with positive energy, at least in the case that the AF interaction 
decays faster then r -4 at large distances. Our argument is substantiated by 
explicit computations in the simple case that the AF interaction depends on 
the Manhattan (L 1 ) distance between sites and decays as r~ p , p > 2, at large 
distances. 

The rest of the paper is organized as follows. In Section 2 we introduce the 
notions of line and corner energies and present our argument explaining why 
stripes are favored as compared to checkerboard when the AF interaction decays 
at infinity faster than r~ 4 and the scale becomes very large compared to the 
lattice spacing. In Section 3 we present detailed analytical computations of the 
stripe and checkerboard energies in cases where the AF interaction depends on the 
Manhattan distance between sites and decays either exponentially or as a power 
law r~ p , p > 2. In Appendix A, we rigorously compute the critical strength 
J c of the FM interaction separating a FM from a non FM phase, when the AF 
interaction decays at infinity faster than r~ 3 . In Appendix B, we prove that power 
law interactions depending on the Euclidean distance between sites are reflection 
positive. This implies that if the GS consists of stripes it will be periodic. Finally, 
in Appendix C we discuss in some more detail the zero temperature phase diagram 
of the model when the AF potential is an exponential Kac interaction: in this case, 
we have evidence for a transition from checkerboard to stripes as J is increased 
from zero to J c . The conjecture is verified by rigorous upper and lower bounds 
on the GS energy. 

2. LINES AND CORNERS 

In this section, we show that the formation of stripes of mesoscopic size in 
d — 2 is essentially independent of the nature of the AF interaction in Eq.(l.l), 
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as long as it is long range and falls off faster than |x|~ 4 , i.e., < i>(x) < A'|x|~ 4 ~ <5 
for some constants K, 5 > 0. According to this argument, the occurrence of 
stripes is related to the sign and the relative sizes of line and corner energies, 
which we now define. 

Consider a system in a square box of side length L with half the spins up and 
half down, separated by a vertical line, called an anti-phase boundary. When the 
falloff of v is faster than |x| -3 the energy divided by L will have a nice limit as 
L — > oo, which is defined to be the line energy r: 

r = -2 1imL- 1 £ £ 0(x-y). (2.1) 

-L/2<a;i<0 l<yi<L/2 
l<x 2 <L l<J/ 2 <i 

The energy per unit length r has the interpretation of surface tension of an infinite 
straight line, and is linear in J, i.e., r = 2(J — J c ) for a suitable positive constant 

Jc- 

At J = J c , the surface tension of an infinite straight line vanishes and there 
is coexistence of the FM ground state with the ground state corresponding to a 
single isolated anti-phase boundary. It is intuitive that for all J > J c , the ground 
state is ferromagnetic, since the energies of ferromagnetic contours (or, at least, 
of straight FM contours) is positive. See Appendix A for a proof of stability of 
the FM state against arbitrary contours. For J < J c the GS is certainly not FM, 
because the system reduces its energy by producing anti-phase boundaries. 

Next, we define a corner energy, k, by first taking two crossed, vertical 
and horizontal, anti-phase boundaries in the box of size L. The energy of this 
configuration is, to first approximation, 2tL. The difference between this energy 
and 2tL has a limit as L — > oo whenever the falloff of v is faster than |x| -4 . This 
difference is the corner energy k, and is given by the formula 

K = 4 E 5>( x -y)+ 4 E 5>( x -y)' ( 2 - 2 ) 

xeQi yeQ3 xeQ2 y€Q4 

where Qi, Q2, Q3, Qa are the first, second, third and fourth quadrant in Z 2 , respec- 
tively. Note that k does not depend on the nearest neighbor interaction energy 
and is therefore positive for the Hamiltonian in Eq.(l.l). 

We now observe that if the GS is made up of rectangles then it necessarily 
consists of a mixture of horizontal and vertical lines, and hence has corners where 
these lines intersect. To lower the energy one can replace the horizontal lines 
by the same number of vetical lines, thereby eliminating the corners. While the 
increased density of vertical lines increases the energy, the saving on the corners 
more than makes up for it when the scale is large enough and J < J c . In 
fact, consider a configuration of sparse straight lines, all at a mutual distance 
larger than R » 1. The interaction energy of any vertical (resp. horizontal) 
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line in a square box of side L with all the other vertical (resp. horizontal) lines 
is positive and smaller than (const.) L Rr l ~ 5 , which follows from the fact that 
< v (x) < i^|x| _4_<5 . Similarly, the interaction energy of any corner with all 
the other corners is negative and smaller in absolute value than (const.) R~ 5 . 
Therefore, the total energy E\ of a configuration of widely separated straight 
lines in a square box A C I? of side L has the form 

E a =(t + OiR- 1 ' 6 )) (M 1 + M 2 )L + (k - 0{R- 5 )) M X M 2 , (2.3) 

where M\ (M 2 ) is the total number of horizontal (vertical) lines. Eq.(2.3) shows 
that, for given M = M x + M 2 of order L, it is energetically favorable to have 
M 1 M 2 = o(L 2 ). In fact, if the number of corners had a finite density, than we 
could decrease the energy by rotating all the vertical (horizontal) lines by 90°, 
making them horizontal (vertical) and placing them half-way between the existing 
horizontal (vertical) lines. After this flipping, the final configuration would have 
an energy equal to (r + 0(R~ 1 ~ 5 )) (Mi + M 2 )L, which is strictly smaller than the 
one of the initial configuration. 

In this sense, corners play the role of elementary excitations, with a positive 
energy cost, which can be eliminated by rotating straight lines by 90°. A similar 
analysis shows that also the "half corners" produced each time that a non-straight 
anti-phase boundary has a 90° turn have a finite positive cost. We are, however, 
not able to exclude the presence of more complicated "excitations" in the GS. 

Regarding the condition on the large distance decay of the AF interactions, 
we do not think it is sharp. However, in the general case, the balance between 
the corner and line energies is much more subtle. In fact, if the decay of the AF 
potential is ~ r~ p , 2 < p < 4, then the corner energy is formally infinite; however, 
corner-corner interactions have an oscillatory sign and such oscillations make the 
effective energy of each corner finite and approximately proportional to R A ~ P if 
2 < p < 4, where R is the distance to the neighboring corner. It is straightforward 
to check that if the corners have a finite density, then their contribution to the 
specific GS energy is comparable to the line-line interaction and of the order of 
R 2 ~ p , where R is the typical separation between lines. Therefore, by rotating the 
vertical lines by 90°, we gain the corner energies and lose some line-line interaction 
energy, both of the order (const. )R 2 ~ P L 2 ; to decide whether the saving makes up 
for the loss, we need to compute the constant prefactors. This will be done 
analytically in the next section, in the special case of AF interactions depending 
on the Manhattan distance between sites. The computation shows that when we 
rotate the vertical lines by 90° and eliminate the corners, the saving overcomes 
the loss for all p > 2. It remains to be seen whether this saving is an accident of 
the specific model considered below or whether there is a general physical reason 
behind the result. 
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We note that in the special case that the AF interaction is reflection positive 
[8], given that the configurations entering the GS are all straight vertical (hori- 
zontal) lines, then they have to be periodically arranged. This follows from the 
analysis in [10-13]. 

3. COMPARISON OF THE STRIPE AND CHECKERBOARD 

ENERGIES 

In this section we perform explicit analytic computations of the energies of 
the stripe and checkerboard states, for different choices of the fall off of the long 
range AF potential. We focus on the (analytically) simple case of interactions 
depending on the Manhattan (L x ) distance ||x||i := |#i| + \x 2 \ between sites. Our 
calculations complement and simplify those in [21]. 

Let us consider Eq.(l.l) with d = 2, e = 1 and 

POO 

u(x) = / dap(a)e~ aMl , (3.1) 
Jo 

with /i(a) a positive measure. We will be particularly concerned with two cases: 

1. Exponential interactions, v(x) = ^e -7 ^ 1 , corresponding to the choice 
p(a) = 7 2 5(a — 7) in Eq.(3.1); 

2. Power law interactions, i>(x) = ||x||^ p , with p > 2, corresponding to the 
choice p(a) = a p - l /T(p) in Eq.(3.1). 

As mentioned above, the choice Eq.(3.1) is made to simplify the computations; 
choosing euclidean rather than Manhattan distance should not make a difference 
from the physical point of view. Let us remark that the potential in Eq.(3.1) is 
reflection positive [8] and so is the (more usual) power law potential v(x) = |x|~ p , 
with |x| = \/x\ + x\ the Euclidean distance (see Appendix B). The property 
of reflection positivity is not explicitly used in the computations below but, as 
observed at the end of previous section, it implies that if the GS consists of stripes, 
then these must be regularly spaced, see [10-13]. 

Let syXx) be the ID profile of period 2h, obtained by extending periodically 
over Z the function / : (—h, h] — > M such that f(x) = sign(x — 1/2) for x = —h + 
1, . . . , h. Let e c (h) be the specific energy of the checkerboard configuration, cr x = 
Sh(xi)Sh(x2), let e s (h) be the specific energy of the striped configuration o" x = 
Sh(xi). We start by computing the specific energy e{h\, I12) of the "rectangular" 
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configuration Sh 1 (xi)sf l2 (x 2 )- We have: 



e(h 1 ,h 2 )= (3.2) 
2J 2J 

~7h 



+ - da^a) E e-l'^le-^-^lx^ ^ <r y ) , 

^ 2 ^ 2 ^ i££ hl ^ 

l<X2<h2 

where x (condition) is = 1 if the condition is satisfied, and = otherwise. After 
some straightforward algebra, 

„ , x 2J 2J f°° da , . 
e(h 1 ,h 2 ) = — + — + 2 — fi(a)- (3.3) 

«1 ^2 Jo a 

tanh(a/ii/2) tanh(a/i 2 /2) tanh(a/i 1 /2) tanh(a/i 2 /2) 
ahi/2 " ah 2 /2 ahi/2 ah 2 /2 

where 

(q/2) 3 cosh(q/2) (a/2) 4 

sinh 3 (a/2) ' a sinh 4 (a/2)- 1 J 

Note that, for small a, A Q ~ 1 - (l/15)(«/2) 4 and 5 a ~ 1 - (2/3)(a/2) 2 , which 
will be useful in the following. 

Using Eq.(3.3), we see that the energy of a striped configuration of period h 
is equal to 

e.ih/2) = £ + 2 f° ^ „(„) [ " 2 A J^M] , (3.5) 
a J cr L ah/2 J 

while the one of a checkerboard configuration of period 2h is 



4J „ f°° da 



h In OL 



tanh(a/i/2) tanh 2 (a/i/ 2 ) 



e c {h) = — + 2\ — /i(«) -2A a \ J ' + B a , _ ^ ' . (3.6) 



a/i/2 (a/i/2) 5 



It is interesting to note that the various terms in Eqs.(3.3)-(3.5)-(3.6) have a 
clear interpretation in terms of the notions of "line energy" and "corner energy" , 
introduced in Section 2. In fact, looking at Eq.(3.3), the terms proportional 
to J correspond to the FM surface tension energy; the integral terms with the 
integrand proportional to A a correspond to the AF line energy (including both the 
negative AF surface tension and the repulsive line-line interactions); the integral 
term with the integrand proportional to B a corresponds to the AF corner energy 
(including both the positive corner self-energy and the attractive corner-corner 
interactions). The analogous terms in Eqs.(3.5)-(3.6) have a similar intepretation; 
note that e c (h) includes a positive contribution from the corner energy, which does 
not appear in e s (h/2), while the contribution from the line energy is smaller than 
the corresponding one in e s {h/2). 

As we discussed above, the goal is to find the balance between these terms 
when the scale of the relevant structures is large compared to the lattice spacing. 
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We will in fact show that when /j>1, then e c (h) > e s (h/2), which is equivalent 
to 



1 r 

2 Jo 



da , . tanh 2 (a/i/2) f°° da . . tanh(a/i/2) - tanh(a/i/4) 
n{a)B a . ,V ; > / — 



(3.7) 

implying that the GS is striped. This will be proved below by treating separately 
the cases of exponential interactions and of power law interactions, with p > 4, 
p = 4, 3 < p < 4, p = 3, 2 < p < 3 (that are listed here in the order of increasing 
difficulty) . 

Remark. Even though Eq.(3.7) does not involve the parameter J, the 
condition that the scale h of the GS structures is large compared to the lattice 
spacing is satisfied only if J is chosen properly. More precisely, as discussed in 
Section 2 (see also Appendix A), if the AF interaction decays faster than r~ 3 , 
then there exists a finite J c such that the homogenous FM state is the GS for all 
J > J c ; in this case, the condition that h ^> 1 is valid in the range J < J c . On 
the contrary, if the decay of the AF interaction is equal to r~ 3 or slower, then 
the condition h 3> 1 is verified for all J 3> 1. The results below are relevant for 
J belonging to these ranges. 



A. Exponential interactions and power laws with p > 4 

If the AF interaction decays exponentially or as a power law with p > 4, 
then we already know from the analysis in Section 2 that e c (h) > e s (h/2) for all 
h ^> 1. For completeness, let us check this analytically, using Eq.(3.7). In the 
case of exponential interactions, the condition reduces to 

1 tanh 2 ( 7 V2) tanh( 7 V2) - tanh( 7 fe/4) 

2 7 ( 7 /i/2) 2 7 lh/2 1 ' ' 

which is obviously satisfied for h large, simply because the l.h.s. goes to zero as 
h~ 2 , while the r.h.s. goes to zero exponentially fast in h. In the case of power 
law interactions with p > 4, the l.h.s. of Eq.(3.7) can be rewritten as 

while the r.h.s. is 

f°° p _ 3 tanh(a h/2) - tanh(a/i/ 4 ) 
/ ua a A a — 
Jo ah/2 

= J- [°° da a^ tanh(a/2)-tanh(q/4) J_ 
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where in estimating the error term of order Yr p ~ 2 we used the fact that | A a — 1| < 
Ca 4 , for a suitable constant C. Therefore, Eq.(3.7) is valid, simply because 
hr 2 > h~ p+2 , Vp > 4, for h large. 

B. The case p = 4 

This case is very similar to the previous one. In fact, the r.h.s. can be rewritten 
and estimated exactly as in Eq.(3.10), with p = 4; in particular, it is ~ h~ 2 . The 
l.h.s. can be rewritten as 

If 00 , „ tanh 2 (a/i/2) 2 f 1 da , 2/ , /n , ^. 1 



2 



/•°° , n tanh^a/i/2) 2 f 1 da , 2/ , /n , ^. 1 , 
/ rfaa£ a \ - = T2 / tanh (a/i/2) + 0(— ) 

Jo (ah/2) 2 h 2 J 1/h a h 2 



i/h 

= 2^ + 0(l). ,,U) 
Therefore, Eq.(3.7) is valid, simply because h~ 2 \ogh 3> /i~ 2 , for ft, large. 



C. The case p < 4 

If p < 4 the proof of Eq.(3.7) is slightly more subtle, because both sides of the 
inequality scale in the same way as h — > oo. In fact, the l.h.s. can be rewritten 

as 



a- 3 ^^ + 0(1] 



a 2 /i 2 ' 

(3.12) 

while the r.h.s. reads 



Jo 

/2\p-2 r 00 tanha - tanh(a/2) 1 , /0 1 oX 

= (s) I &Q ^ t0 W' (3l3) 



p-3 a tanh(a/i/2) — tanh(a/i/4) 
cto A a - — 

cm/2 



Therefore, both sides of Eq.(3.7) scale as ~ h 2 p as h — > oo. The inequality is 
asymptotically valid if and only if the following condition is true: 



1 r 

2 Jo 



„ , tanh 2 a f°° , „ , tanha - tanh(a/2) 

daa p 6 5 — > / daa p s ^ . (3.14) 

« 2 Jo « 



This inequality can be checked numerically in the different ranges 3 < p < 4, 
p = 3 and 2 < p < 3. In fact, if p = 3, Eq.(3.14) is equivalent to 



f°° tanh 2 a 

/ cfa — = 0.85256... > log 2 = 0.69315... (3.15) 

2 Jo « 
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If 3 < p < 4, Eq.(3.14) is equivalent to 

1 roc roc 

- da a p ~ 5 t&nh 2 a > (2 p - 3 - 1) da a p " 4 (l - tanh a) . (3.16) 

2 Jo Jo 

In the limiting case p — > 3 + , Eq.(3.16) is equivalent to Eq.(3.15), as it should. 
In the limit p — > 4~, Eq.(3.16) is obviously valid (because the r.h.s. tends to a 
constant, while the l.h.s. diverges to +oo). The validity of Eq.(3.16) for all values 
of p in the interval (3,4) can be checked numerically, see Fig.l. 

6 - 

5 - 

4 - 

3 1 

2 - 




3.0 3.2 3.4 3.6 3.8 4.0 

FIG. 1: A plot of the difference between the right and left hand sides of (3.16) vs p, 
which proves that e s (h/2) < e c (h) for all 3 < p < 4, and h large enough. 

Finally, if 2 < p < 3, Eq.(3.14) is equivalent to 

^ roc poo 

- da a p ~ 5 tanh 2 a > (1 - 2 P ~ 3 ) / da a p ~ 4 tanh a . (3.17) 
2 Jo Jo 

In the limit p — > 3~, condition Eq.(3.17) reduces to Eq.(3.15), as it should. In 
the limit p — > 2 + , condition Eq.(3.17) reduces to 

f°° r /tanha tanh 2 a\ .„ . 

log2> / da( — ^ 3.18 

Jo v a a ' 

^ log2 + - = 1.193147... > / rfa tanl ! - = 1.154785... 
2 J a 2 

The validity of Eq.(3.17) for all values of p in the interval (2,3) can be checked 
numerically, see Fig. 2. 

This concludes the proof that e c (h) > e s (h/2) whenever h is large, for all 
power law decays with exponent p > 2 and for exponential interactions. An 
immediate consequence of this analysis is the following: let e*(J) = min he ^ e s (h) 
and e*(J) = min heN e c (h) be the optimal stripe and checkerboard energies at 
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2.0 2.2 2.4 2.6 2.8 3.0 

FIG. 2: A plot of the difference between the right and left hand sides of (3.17) vs p, 
which proves that e s (h/2) < e c (h) for all 2 < p < 3, and h large enough. 

a given J; then, if the AF interaction is either exponential or power law with 
p > 3, we have e*(J) < e*(J) for all J c — J positive and small enough; if the 
AF interaction is power law with 2 < p < 3, then e*(J) < e*(J) for all J large 
enough. 

In conclusion, we showed for a 2D spin model with competing short range 
(nearest neighbor) FM and long range AF interactions that stripes are favored 
with respect to checkerboard when the GS structures are large compared to the 
range of the FM interaction. If the AF interaction decays faster than r -4 , the 
emergence of stripes close to the transition to the FM phase can be understood 
on the basis of a comparison between the sign and relative sizes of the corner 
and line energies, which is independent of the details of the AF interaction. If 
the decay at infinity of the AF interaction is slower, than the balance between 
corner and line energies is more subtle, and the understanding of why stripes are 
favored relies on explicit computations of the stripe and checkerboard energies, 
which have been performed here in the simple case that the AF depends on 
the Manhattan distance between sites. We believe that future progress on the 
problem will come from a deeper understanding of the reason that interactions 
that fall off slower than r~ 4 always seem to favor stripes. 
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Appendix A: A rigorous computation of J c 

Let us assume that the long range AF interaction decays at infinity faster 
than r~ 3 , and let r = 2(J — J c ) be the line energy, as defined in Section 2. Here 
we want to prove that for all J > J c , the homogeneous FM state is a GS of 
Eq.(l.l) (and is the unique GS for J > J c ). As already remarked in Section 2, 
for J < J c the homogenoeus state is not a GS, simply because the state with a 
single straight anti-phase boundary has negative energy. If J > J c we want to 
get a lower bound on the energy of an arbitrary state, which is positive, unless 
the state is homogeneous. 

We proceed in a way similar to the proof of Theorem 3 of [10]. We need to 
introduce some definitions; in particular via the basic Peierls construction we 
introduce the definitions of contours and droplets. Given any spin configuration 
cr A on the squared periodic box A, we define A to be the set of sites at which 
<Tj = —1, i.e., A = {i G A : (Xj = —1}. We draw around each i 6 A the 4 
sides of the unit square centered at i and suppress the sides that occur twice: we 
obtain in this way a closed polygon T(A) which can be thought as the boundary 
of A. Each side of T(A) separates a point i G A from a point j G" A. At every 
vertex of T(A) n A*, with A* the dual lattice of A, there can be either 2 or 4 sides 
meeting. In the case of 4 sides, we deform slightly the polygon, "chopping off" 
the edge from the squares containing a — spin. When this is done T(A) splits into 
disconnected polygons r 1; . . . ,T r which are called contours. Note that, because 
of the choice of periodic boundary conditions, all contours are closed but can 
possibly wind around the box A. The definition of contours naturally induces a 
notion of connectedness for the spins in A: given i, j G A we shall say that i 
and j are connected iff there exists a sequence (i = io,ii, ■ ■ ■ ,i n — j) such that 
im, im+i, m — 0, . . . ,n — 1, are nearest neighbors and none of the bonds (i m , i m +i) 
crosses T(A). The maximal connected components Si of A will be called droplets 
and the set of droplets of A will be denoted by T>(A) = {Si, . . . , S s }. Note that 
the boundaries T(Si) of the droplets Si G T>(A) are all distinct subsets of T(A) 
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with the property: U s l=1 T(5i) = T(A). 

Given the definitions above, let us rewrite the energy E^(a A ) of q_ h in a box 
AcZ 2 with periodic boundary conditions as 

E A (a A )=2J l r |- E ( A ' X ) 

rer(A) <5ex>(A) 

where E dip (5) := 2e X] x e5 S y eA c t '( x ~~ 30' w hich can be bounded from above as 

E dip (5) = 2eJ2 ^( n ) EE^-^^ 

nez 2 xe<5 yeA c 

< 2eJ2v(n)J2 E X(x-y = n). (A.2) 

nez 2 xe<5 y ez 2 \<5 

Now, the number of ways in which n = (ni, n 2 ) may occur as the difference x — y 
or y — x with x G 5 and y ^ 5 is at most Xlrer(<5) Ya=i |r|i|^i|, where |r|j is the 
number of faces in T orthogonal to the i-ih coordinate direction. Therefore, 

2 

E dip (5) < ^ E E v ^ E i r i*H = 2e E i r i E nit; ( n ) = 2Jc E i r i • 

rer(5)nez 2 i=i rer(<5) n ez 2 : rer(<5) 

ni>0 

(A.3) 

Plugging this back into Eq.(A.l) gives 

E A {a K )>2{J-J c ) E l r l> ( A - 4 ) 
rer(A) 

which readily implies that the uniformly magnetized state is a GS for all J > J c 
and that it is the only GS for J > J c . 

Appendix B: Reflection positivity of power law interactions 

In this Appendix we prove that v (x) = |x|~ p , with p > and |x| = |x| 2 = 
\fx\ + x\ the usual Euclidean distance, is a reflection positive (RP) potential, 
which may be a useful remark for a possible future proof of the periodicity of the 
GS of Eq.(l.l) with v (x) = |x| _p . We recall that v is RP if, for all compactly 
supported functions / : Z 2 — > C, 

E /x/ y ^i+yi-M2-!te) >o. (B.i) 

X!,y!>l 

By Schur's product theorem, pointwise products of RP potentials are reflection 
positive. Therefore, in order to prove that |x|~ p is RP for all p > 0, it is enough 
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to show that |x| 1 and |x| A , with < A < 1, are separately RP. If v (x) = |x| 1 
and X\ > 1, 



v(x 1 ,x 2 ) = , = = 7TT / dk I d P d 1 



y/xj + x\ 2tt 2 k 2 +p 2 + q 2 

= - / dp dq e-ziV^W , (B.2) 

^ J-oo J-oo VP 2 + q 2 

from which (B.l) readily follows. If v (x) = |x|" A , with < A < 1, then (B.l) 
follows if we prove the stronger result 

dx x f dy x [ + dx 2 ^ dy 2 P ,^ P ^} > , (B.3) 



i x -yi 



if p(x) is a smooth compactly supported real function, with support contained in 
R 2 \ {xi = 0}, and such that p(—Xi,x 2 ) = p(xi,x 2 ). Using the Fourier transform 
of |x|~ A , see e.g. [18] Thm. 5.9, and proceeding as in [7], we can rewrite the l.h.s. 
of (B.3) as 

I Y(l — -) f f e iki(x!+yi) e ik 2 (x2-y2) 

^^W~ U****"™ WTW^ p(y) ' (B ' 4) 

We observe that for fixed x 1 +y 1 > and k 2 , the function e ikl ( Xl+yi \k 2 + k?,y 1+x/2 
is analytic in k\ in the upper half plane with the cut {ir : r > \k 2 \} removed. 
Deforming the contour of integration in dk\ to this cut and calculating the jump 
of the argument across it we obtain 

/+oo gifci(xi+j/i) /-oo g -T(xi+j/i) 

.„ ^' PF = 2 (7F(1 " V2)) i 2 , dT (r»-ifci)i-V2 • ( B - 5 ) 

Plugging this back into (B.4) we find 

- 2 — F J^ sm („(1 - A/2)) / dk2 f* . (B.6) 

JR J\kn\ 



I- 



-*1.W>° |x-y| A 7T r(f) 7 J K i| fc2 | 

which is clearly nonnegative. This concludes the proof that |x|~ p is reflection 
positive for all p > 0. 



Appendix C: Kac interactions 



In this appendix we add some comments about the possible structure of the 
GS of Eq.(l.l) in the case that v is a 2D Kac potential, i.e., i>(x) = 7 2 t> (7x), 
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with 7 a small parameter. These may be relevant for the understanding of the 
"froth problem", addressed by Lebowitz and Penrose in [17]. To be definite and 
make things simple, we restrict to the case of exponential interactions depending 
on the Manhattan distance: v (x) = 7 2 e~ 7 'l x 'l 1 . In this case, J c = 2'j~ 1 A 1 and, if 
J > J c , the GS is the homogeneous FM state. 

From the computations in Section 3, we already know that, as J — > J~ , 
the stripe state is energetically favored as compared to the checkerboard state. 
If e*(J) = min/jgpj e s (h) and e*(J) = mm h€ ^e c (h) are the optimal stripe and 
checkerboard energies, an explicit computation shows that, if < £ := 7(J C — 
J) « 1, 

correspondingly, the scales h* and h* of the optimal stripe and checkerboard 
configurations turn out to be: 

?: == = 2 -f L + 0(e'^) , q* s :=^f = \\ log£| + 0(\og | log^) , (C.2) 

Using methods similar to those in the proof of Theorem 3 of [10], it is easy to 
prove that for £ small the scaling of e*(J) is the optimal one; i.e., the absolute 
ground state energy per site eo admits a lower bound of the form eo > — (const. )£• 

It is interesting that the model with exponential Kac interactions also displays 
a phase where mesoscopic checkerboard are energetically favored with respect to 
stripes. In fact, note that the periods of the optimal checkerboard and striped 
states are given by 

7 J = (2A, - 25 7 ^^) [tanhg c * - q* c (l - tanh 2 q* c )} , (C.3) 
7 J = 2A 7 [tanhg;-g;(l-tanh 2 g;)] , (C.4) 

from which we immediately recognize that, if 7 <C 1 and 1 < J <C 7 _1 , then 
h* and h* are both <C 7 _1 ; therefore, the solution to these equations can be 
determined by expanding their r.h.s. in Taylor series in q and solving to dominant 
order, which leads to: 

2£ = + 0((7 ,p) , ^ -2 + + OfrJf*) ■ 

f = (^) V3 + 0(7 J) , < = -2 + 2(^) 2 ' 3 + 0(( 7 J)V.) . 

Therefore, in this regime the specific energy e* of the optimal checkerboard config- 
uration is smaller than the specific energy e* of the optimal striped configuration. 
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This suggests that for any fixed J and 7 small enough the ground states of the 
considered model display periodic checkerboard order, a conjecture supported by 
the fact that the absolute ground state energy per site admits a lower bound of 
the form eo > —2 + (const.) (7 J) 4 / 5 , which has the right scaling, see below for a 
proof. 

In conclusion, if the AF interaction is exponential with a Kac-like scaling, we 
expect that as J is increased from to J c , the GS should display a transition from 
checkerboard to stripes. On the basis of the previous computations, we expect 
the trabnsition to take place at values of 7J of order one, see Fig.3. 




FIG. 3: A plot of the optimal checkerboard energy e* (solid line) and of the optimal 
striped energy e* (dashed line) vs J := 7 J for exponential Kac interactions v(x) = 
7 2 e _7 " x " 1 at 7 = 0.4. The plot shows a transition from a case where e* < e* (for small 
values of J) to a case where e* > e* (for larger values of J). In the limit 7 <C 1, the 
transition is expected to occur for J of the order 7 _1 . 

Remark. The scaling of the checkerboard energy as well as the very existence 
of a checkerboard phase may depend on the specific choice of the Kac potential. 
In particular, it may depend on the reflection positivity property of the Kac 
potential (note that the considered exponential interaction is reflection positive); 
if vq is smoother at the origin (e.g., t>o(x) = e~' x ' 2 ), the checkerboard phase 
may disappear or, at least, be characterized by a completely different scaling 
behavior. The reason for this is already apparent in a ID toy model: consider 
model Eq.(l.l) in d = 1 with v(x) = 7^0(72;) and Vq either of the form Vq(x) = 
e -N Qr Vq ( x ^ — e~ x2 ; if one optimizes the energy of a configuration consisting of 
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blocks of uniformly magnetized spins of size h and alternating sign, the optimal 
size turns out to be of the order 7~ 2 / 3 in the exponential case and 7~ 1 /l log7| in 
the gaussian case. This can be seen as follows: the scale of the optimal periodic 
structure can be found by balancing the energy contributions from the FM and AF 
interactions; while the first is 2J/h, the second is of the order of Vo(l/~fh), with vq 
the Fourier transform of Vq] the latter depends on the smoothness properties of Vq 
and, more specifically, it behaves like vo(k) ~ k~ 2 or ~ e -(const)/c 2 i ar g e j n ^he 
cases of v exponential or gaussian, respectively. Minimization of 2J/h + v (l/ r yh) 
over h gives the optimal size of the structures. 

The fact that the nature of the checkerboard structure depends on the reflec- 
tion positivity properties of the Kac potential is consistent with the fact that the 
proof of the lower bound on the energy in the Kac regime heavily uses reflection 
positivity, see next subsection. 



1. Lower bound on the energy: Kac regime 



Let us assume that 1 < J <C 7 -1 : in this case we want to prove that 
eo > —2 + (const.) (7 J) 4 / 5 , which asymptotically matches the upper bound e < e* 
and supports the conjecture that, in this regime, the ground state has checker- 
board order. Let E A (a A ) be the energy of the spin configuration g_ A in the peri- 
odic squared box A. Let us consider a partition of A into squares Qi of side I: 

IAI /^ 2 

A = U 

=1 ; given q_ A and Qi we shall denote by <7q. the restriction of the spin 
configuration a A to the square Qi. Let f^(x) = 7 2 Xl ne ^ 2 e 7 ' x+nL ' 1 an d let us 
rewrite 
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Ea(<L a ) = -2 ta ^ /2) |A| + E*(a A ) + E%a A ) , (C.5) 
where E$(a A ) = ±£ 

xjeA^^ — y)°x0y is the antiferromagnetic energy associ- 
ated to cr A , while Ej(q_ A ) = 2J£ xeA Y^=i x(°x 7^ °x+ej is the surface tension 
energy of a A in the box A with periodic boundary conditions. If we drop the 
surface tension energy across the boundaries of the squares Q iy we get a lower 
bound on the energy of the form: 

(>y/2V |A|/£2 ~ 

Ea(£a) " " 2 tan^(7/2) |A| + ^ + ^ ' (C ' 6) 

where Ej^aqJ is the surface tension energy of the spin configuration <7q. in the 
box Qi with open boundary conditions. If Wj := £~ 2 £xeQi(— Q») x ' ^ ne surface 
tension energy can be further bounded from below by: 

Ef(g_ Qi ) > 2J^min{l, 2^/2(1- \m~[j} . (C.7) 
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Moreover, using reflection positivity [8], the antiferromagnetic energy can be 
bounded from below as 

\k\/e. 2 
i=i 

where e 7 (q_Q.) is the specific energy of the infinite volume configuration obtained 
by repeatedly reflecting a_ Q . (with "antiferromagentic reflections" ) across the sides 
of Qi and its images. More explicitly, 

2 ^ A 1 



^)= ? e i- P iv(i - e-) 2 n (1 _ e _ 7)2 + 2e _ 7(1 _ C0SPi) , 

P=f (ni,n 2 ) t=l V ; V Fi; 

rii=l,3,5,...,2£-l 

(C.9) 

where <r p := XlxeQi °xe~ lkx - Using the fact that, for all e > 0, 

ISpl 2 > 7^ ^ r - V(l - \ mi \f , (CIO) 

(1 — cospi)(l — cosp 2 ) £ 



we get 

e 7 (^) > (1 - £)e 7 (l Qi ) - (const.) 1(1 - K|) 2 ( 7 Q 4 , (C.ll) 

where e 7 (l Qi ) is the antiferromagnetic energy per site of the checkerboard config- 
uration with tiles of side £. Note that e^lq.) scales as (const.) (7^) 4 in the regime 
under consideration and for £ ^> 1; moreover, it can be bounded from below by 
C^iy for a suitable constant C. Optimizing over e we get 

e 7 (^) > e 7 (l Qi ) - c(l - K|)( 7 ^) 4 , (C.12) 

for a suitable constant c. Combining all the previous bounds we find that 
Eh{g_ A ) + 2 (7/2) 2 tanh~ 2 (7/2)|A| can be bounded from below by 

f E It min{1 ' 2v/2(1 " Kl)} + P " c(1 _ Kl) ] (y)4 / • (c - 13) 

i=l 

Optimizing over m; and £ leads to £ = (const.) (J7"" 4 ) 1 / 5 and 

eo > -2 ta ^ 2 /2) |A| + (const.)( 7 J) 4 / 5 , (C.14) 

as desired. The proof of (C.14) can be easily adapted to higher dimensions and 
to cases where the ferromagnetic interaction has finite range rather than being 
nearest neighbor. On the contrary, the assumption of RP was used in a crucial way 
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and it is likely that in the presence of more general long-ranged antiferromagnetic 
interactions the ground state energy scales differently with 7. 
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